Spin-triplet 'f-wave' pairing proposed for an organic superconductor (TMTSF) 2 PF 6 
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By examining how the spin- and/or charge-fluctuation exchange can contribute to pairing insta- 
bilities, we propose that a spin-triplet f-wave-like pairing with a d-vector perpendicular to the &-axis 
may be realized in (TMTSF^PFg due to (i) a quasi-one-dimensional Fermi surface, (ii) a coexis- 
tence of 2kF charge fluctuations and spin fluctuations, and (iii) an anisotropy in spin fluctuations. 
Fluctuation-exchange study for the Hubbard model confirms the point (i), while a phenomelogical 
analysis is given for (ii) and (iii). The proposed pairing is consistent with various experiments. 
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Spin-triplet pairing is conceptually fascinating, but 
there seem to be few examples. Recently an organic 
superconductor, (TMTSF) 2 PF 6 , has attracted much at- 
tention since a triplet pairing is suggested from an ob- 
servation of large H C 2 [jj as well as from a Knight shift 
experiment || , while the absence of Hebel-Slichter peak 
and the power-law decay of T-T" 1 below T c || suggest an 
anisotropic pairing with nodes in the gap. 

If triplet pairing is indeed realized in (TMTSF) 2 PF 6 , 
its mechanism is a challenging theoretical puzzle: in the 
pressure-temperature phase diagram for this material the 
superconductivity lies right next to the 2kF spin density 
wave (SDW), so that if one seeks an electronic origin, 
a spin-singlet d-wave-like pairing mediated by spin fluc- 
tuations is most naturally expected as proposed by sev- 
eral authors. |^-|| If on the other hand one assumes a 
phonon-mediated attractive interaction, triplet pairing, 
with nodes in the gap in general, would seem less favor- 
able compared to singlet s-wave pairing without nodes. 
Recently, Kohmoto and Sato Q have proposed that this 
difficulty in the phonon mechanism may be circumvented 
for TMTSF compounds, where they show that the quasi- 
one dimensionality of the Fermi surface, along with the 
presence of spin fluctuations, makes a triplet p-wave pair- 
ing without nodes on the Fermi surface jjj dominate over 
s-wave pairing. However, it is not clear if such a nodeless 
gap can be reconciled with the absence of Hebel-Slichtcr 
peak and the power-law Tf 1 in (TMTSF) 2 PF 6 . 

If we turn to another prominent candidate for triplet 
superconductivity accompanied by SDW fluctuations, ||] 
Sr2Ru04, Takimoto recently proposed that charge fluc- 
tuations (or more precisely orbital fluctuations) should 
arise from repulsions between degenerate 4d orbitals, and 
that the coexistence of spin and charge fluctuations may 
lead to a triplet pairing. J^] This makes us recall an ex- 
perimental fact that a Ikp charge density wave (CDW) 
actually coexists with the SDW in (TMTSF) 2 PF 6 as sug- 
gested from X-ray diffuse scattering. [ |Ic| I n another the- 
ory for Sr 2 Ru04, Sato and Kohmoto, and indepen- 
dently Kuwabara and Ogata p3], have proposed that 
anisotropy of the spin fluctuations, known to be present 



experimentally, |Tg| ] may give rise to a triplet p-wave pair- 
ing. The anisotropy of spin fluctuations is also present in 
(TMTSF) 2 PF 6 . |1 Moreover, Sr 2 Ru0 4 has two quasi- 
ID Fermi surfaces (although they are weakly hybridized 
to result in two 2D Fermi surfaces), so the ruthenate 
seems to share several features with the TMTSF com- 
pound, although the strong charge fluctuation employed 
in Takimoto's mechanism is yet to be detected experi- 
mentally in Sr 2 RuC>4. 

However, the triplet pairing mechanism of 
(TMTSF) 2 PF 6 cannot be the same with that of Sr 2 Ru0 4 
since d (the <i-vector characterizing the triplet pairing) 
_L z (easy axis of the spins) is experimentally suggested 
in the former, while d || z in the latter. |y| In the 
present paper, we propose that a triplet f-wave-like pair- 
ing with d L z can take place in (TMTSF) 2 PF 6 due to 
a combination of (i) the quasi-one-dimensionality of the 
Fermi surface, (ii) coexistence of 2fcp spin and charge 
fluctuations, and (iii) the anisotropy in the spin fluctua- 
tions. In the first part of the paper, we focus on how the 
quasi-one-dimensionality works favorably for the triplet 
pairing, and perform a fluctuation-exchange (FLEX) pl| 
calculation for the on-site repulsion Hubbard model on 
a lattice for (TMTSF) 2 PF 6 . Then, in the second part, 
we discuss phenomelogically how the triplet pairing can 
become competitive against the singlet when charge fluc- 
tuations coexist with spin fluctuations. We finally point 
out that the anisotropy in the spin fluctuations should 
further favor triplet pairing with d _L z . 

We first consider the on-site Hubbard model, TL = 
X)<i,3)o- *»i c i<T C j<r + UJ2i n iT n ili in tne h° le picture on a 
quasi- ID lattice (\ts\ > \ti\) depicted in Figjl|. There are 
n = 0.5 holes per site. Since sites A and B are inequiv- 
alent for tsi ^ t$2 and ti\ ^ ti% (dimerization of the 
molecules), we adopt the two-band version of the FLEX 
p8| , ^9[ (although we shall see that the dimerization is not 
essential in our argument). 

For later discussions, we first recapitulate the one-band 
version of FLEX in a general fashion, where we proceed: 
(i) Dyson's equation is solved to obtain the renormalized 



1 













— ) 

/ 


13 

/ \ 


12/ 
f B t 






b 


' SI / 

11/ 13 


' S2 / 

















-> a 



FIG. 1. The lattice considered with the hopping integrals 
taken to be t si = -2.8, t S2 = -2.5, t n = +0.2, t I2 = +0.5, 
to = —0.5 in units of a typical energy scale, 100 meV, for 
organics after ref. Jit]] 



Green's function G(k), where k is a shorthand for the 
wave vector k and the Matsubara frequency, ie n , (ii) the 
fluctuation-exchange interaction V^(q), given as, |2yj 



V {1) {q) = \v£{q) + V+-(q) + ^(q), 



(1) 



consists of the contribution from longitudinal (zz) and 
transverse (H — ) spin fluctuations (sp) and that from 
charge fluctuations (ch). For the on-site Hubbard model 
in particular, V zz = V+~{= V sp ) = U 2 X sp and V ch = 
U 2 x ch , where the spin and the charge susceptibilities are 
given as X sp (<l). = X to (<?)/[1 - U^(q)} and X ch (g) = 
X irr (g)/[1 + Ux 11T {q)}, respectively, using the irreducible 
susceptibility x irr (<?) = -]? E fc G(k + q)G{k) (iV:number 
of fc-point meshes), (iii) then brings about the self- 
energy, E(fe) = i J2 q G(k - q)VW(q), which is fed back 
to Dyson's equation, and the self-consistent loop is re- 
peated until convergence is attained. 

T c is the temperature where the eigenvalue A of the fol- 
lowing Eliashbcrg equation for the superconducting order 
parameter 4>(k) reaches unity. 



T 

'N 



^2Mk')\G(k')\ 2 vW(k-k>) (2) 



In the two-band version of FLEX, G, \i and <fi all 
become 2x2 matrices, whose elements are denoted as 
G a p etc with a, j3 —A or B in the site representation, 
which may be converted to the band representation with 
a unitary transform. Since the Fermi surface lies in the 
lower band for quarter filling, we concentrate on Green's 
function and the order parameter in that band, denoted 
as G and <f> s , <j>t, respectively. As for the spin susceptibil- 
ity, we diagonalize the 2x2 matrix x sp and concentrate 
on the larger eigenvalue, denoted as x- To ensure con- 
vergence at low temperatures in the two-band system we 
had to take 64 x 64 fc-points and e n from — (2N C — l)7rT 
to (2N C - 1)ttT with N c up to 8192. 

In Fig.||, we show contour plots of |G(k, i7rfcsT)| 2 (a), 
X(k,0) (b), 4> s (k,iTTk B T) (c), and fo(k,iTrk B T) (d) for 
T = 0.015. The Fermi surface as identified from the 
ridge in |G(k)| 2 is a pair of warped quasi ID pieces. The 
spin susceptibility x(q, 0) has a peak at Q ~ (n, it/2) (or 
(7r/2, 7r/2) in the unfolded Brillouin zone in the absence of 
dimerization), as expected from the nesting vector and in 
agreement with experimental results. \21 22] The singlet 



pairing order parameter is seen to change sign in such a 
way that (i) s (k) = S (— k), and (ii) <^ s (k) has opposite 
signs across the nesting vector Q so that the repulsive 
Vs 2 ^(Q) (eq.(||)) acts as an attractive interaction in the 
gap equation. We call the singlet pairing a 'd-wave' in 

that the sign of s (k) changes like H h— if we rota- 

tionally scan the Fermi surface, which is consistent with 
previous studies 

(2) 

For the triplet pairing, by contrast, V t (Q) is attrac- 



tive (eq.(W) or (Rh with V s j 



V zz ), so that the order 



parameter should have the same sign across Q. This re- 
quirement, along with the condititon for a triplet order 
parameter </>t(k) = — 4>t{— k), can be satisfied by adding 
extra nodal lines along k a ~ and k a ~ it (mod 27t). We 
call this pairing an 'f-wave' in that <p s behaves this time 
like H 1 1 — along the Fermi surface. 



Here, the pairing interaction V^i (q) is given as 
V s {2) (q) = \v™(q) + V+-(q) ±V ch (q) 



for singlet pairing, 

V t { l\q)- 



-^V s z p z (q) - -V ch (q) 



(3) 



(4) 



for triplet pairing with total S z = ±1 (d _L z), and 

V t?M = l V sf(q) - V+~(q) \v ch {q) (5) 

for triplet pairing with S z = (d \\ z). In the on- 
site Hubbard model, V sp 3> Kh is satisfied, so that 
\V t {2) \ ~ (1/3)|^ (2) | holds with V t (2) = V t { l } = V t {2) . 
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FIG. 2. Contour plots of |G(k, ink B T)\ 2 (a), x (q, 0) (b), 
(f> 3 (k,iirk B T) (c), and <j) t {^,mk B T) (d) for n = 0.5, (7 = 8, 
and T = 0.015. 
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FIG. 3. At/A s plotted as a function of temperature for 
the parameter values adopted in Figj^ (solid line), or for 
ti2 = —tiz = 1.2 (dashed line). 

A virtue of the quasi-one dimensionality is that the 
magnitudes, |0 s (k)| and t (k)|, are almost identical 
around the Fermi surface as seen from Figs.^(a) and 
(c),(d). In fact, the 'f-wave' here is to 'd-wave' what 
p-wave is to s-wave in Kohmoto-Sato's picture || in that 
a singlet is converted into a triplet by introducing extra 
nodes that do not affect |</>(k)| on the Fermi surface. In 
such a situation, the difference between A s and At comes 

(2) 

almost entirely from the difference between \ V S (Q) and 
v/ 2 '(Q)| in the Eliashberg equation (||). 

Our result shows that the ratio A t /A s indeed tends to 
1/3 at low temperatures (Fig.||; solid line), which reflects 
the ratio \V t {2) (Q) /V S {2) (Q)\ ~ 1/3 in the on-site Hub- 
bard model. The ratio approaches to 1/3 as the temper- 
ature is lowered, because the ridge in \G\ 2 sharpens so 
that the triplet pairing, with the order parameter van- 
ishing around k a — 0, it, becomes more favorable. We 
can also confirm that quasi-lD is exploited in realizing 
At /A s ~ 1/3 by pushing the system toward 2D with larger 
value of ti2 and \ti3\, for which the ratio X s /Xt deviates 
from 1/3 even at low temperatures (FigJ|; dashed line). 

We have so far seen that the difference in A between 
singlet 'd' and triplet 'f ' can directly reflect the difference 

(2) (2) 

between \V S ^'(Q)\ and \V t y >(Q)\ in a quasi- ID system. 
The 'f-wave' proposed here is an appealing candidate, be- 
cause it can account for both experimentally suggested 
triplet pairing and the nodes in the superconducting gap. 
However, even for a quasi-lD system, 'f is only 1/3 com- 
petitive against 'd' as far as the on-site repulsion Hub- 
bard model is concerned — to make 'f dominate over 'd', 
we need to have \V t {2) (Q)\ > |V^ (2) (Q)|. So at this stage 
we depart from the Hubbard model to argue phenome- 
logically how some factors in the actual (TMTSF) 2 PF 6 
that are not taken into account in the simple Hubbard 
model can indeed make 'f competitive against 'd'. 

An important experimental fact for (TMTSF) 2 PF 6 
that cannot be explained by the on-site Hubbard model 
is that a 2kp CDW actually coexists with the 2kp SDW 
. The coexistence of spin and charge fluctuations can 
favor triplet pairing as pointed out for Sr 2 Ru04 by Taki- 
moto mentioned above. || This can be seen in eqs.(|3|), ([!]), 
where an increase in V^h enhances v/ 2 ^| and suppresses 



\V} '\ (as far as V sp > 3V c h for isotropic spin fluctuations 
assumed for the time being). Now, if we take the coex- 
istence of 2k F SDW and 2k F CDW in (TMTSF) 2 PF 6 
to be V sp (Q) ~ Vch(Q), eqs.(f| and (@) dictate that 
\V t {2) (Q)\ ~ \V s {2 \Q)\, so 'f does indeed compete with 
'd', but the competition is still subtle. 

Is there any mechanism that further favors the triplet 
pairing? Magnetic anisotropy is, in our view, one. It has 
actually been revealed experimentally for (TMTSF) 2 PF6 
that the SDW has an easy axis in the ^-direction, |l4| 
which implies that V s z p z (Q) > V^+~(Q) is satisfied for 
z taken to be || b. In such a situation, the 'f pair- 
ing is more favorable in the S z = ±1 channel since 

(2) (2) 

\V t \ \ > \VX |. The condition for 'P dominating over 
'd' now reads \V^(Q)\ > V S {2) (Q), or 

^ch(Q) > V+-(Q) (6) 

from eqs.(||),(^). This last condition should be satisfied 
in (TMTSF) 2 PF6 because the spins do not order in the 
transverse direction even in the SDW phase, while the 
charges do. We stress that d _L z with z || b is consistent 
p3| with the experimental result: it is when the magnetic 
field is applied parallel to the b-axis that (i)iJ c2 becomes 
largest at low temperatures, ||] and (ii) the Knight shift 
is unchanged across T c 

The mechanism in which the anisotropy of the spin 
fluctuations favors triplet pairing is reminiscent of the one 
proposed in refs. OjTJJjfor Sr 2 Ru04, but a crucial dif- 
ference is that refs. ]ll| , p^ | , which do not consider charge 
fluctuations, conclude a p-wave pairing with S z = in 
agreement with the experimental results suggesting d \\ z 
in Sr 2 RuC>4. [jl5| Let us see how this would occur in the 
present context. If V s z p z > {2V<+~ + V c h), we can see from 

cq.(||) that V^y (Q) becomes repulsive, which will medi- 
ate a triplet pairing having an order parameter with op- 
posite signs across Q. This requirement, along with the 
triplet condition 0t(k) = — <f>t{— k), can be satisfied by 
putting nodes only at k a ~ and k a ~ tt, thereby mak- 
ing the order parameter nodeless on the Fermi surface as 
in refs. Specifically, if V ch /V z p z and V+~/V^ z are 

(2) 

both sufficiently small, the pairing interactions V} (fa- 

(2) (2) 

voring 'd'), V tJ _ ('f'), and VX (p) will all have similar 
magnitudes, so the p-wave pairing, with no nodes on the 
Fermi surface, should dominate over the others. Thus, 
the 'f pairing is not realized unless V^h/Kp 2 1S signifi- 
cant even if eq.(^) is satisfied (see Fig.||). 

In this context, a possibly related problem is the pair- 
ing symmetry in (TMTSF) 2 C104, another candidate for 
a triplet superconductor. For this compound an NMR 
experiment suggests a presence of nodes on the Fermi 
surface, J24[] while a recent thermal conductivity 
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FIG. 4. A phase diagram against the charge/spin axis 
(V c h/V S p Z ) and the spin anisotropy axis (14p~ /V s z p z ). The solid 
line is according to eq.(H), while dashed lines are schematic. 
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to second nearest neighbors. |3(| It would be an inter- 
esting future problem to investigate microscopically the 
singlet-triplet competition and to evaluate T c by taking 
these effects into account. 
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measurement suggests a nodeless gap. pq] If we adopt the 
latter result, the nodeless p-wave pairing should become 
a strong candidate. Then a comparison of the magni- 
tude of the charge fluctuation as well as the direction of 
d between (TMTSF) 2 C10 4 and (TMTSF) 2 PF 6 will be a 
crucial test. 

Having discussed the lower bound of V c b./V^ for 'f 
pairing, how about an upper bound? When V c h/V s z p z ^> 
1, singlet s-wave pairing with </> s (k) ~ constant should 



enter as the dominant pairing. This is because V s 
comes attractive for V c h > (V S p 



(2) 



be- 



- 2V+-), so that (f> s (k) 

(2) (2) 

no longer has to change sign. As V~ s /V t \ tends to unity 
with the increase of charge fluctuations, the 'F, with its 
nodes on the Fermi surface, should thus give way to the 
nodeless s. All the above reasoning is schematically sum- 
marized as a generic phase diagram in Fig.^. 

An additional bonus from the coexistence of strong 
spin/charge fluctuations and the anisotropic spin fluctu- 
ations is that these effects may serve to enhance the tran- 
sition temperature for triplet pairing. Namely, a flaw in 
triplet superconductivity mediated by isotropic spin fluc- 
tuations is that the absolute value of the triplet pairing 

(2) 

interaction \ V t | is only one third of the effective interac- 
tion that determines the normal self-energy as seen 
from eqs.(||) and (@) (or (§) with U/ p z = V+~ » V ch . 
This is in contrast with the case of singlet pairing, where 
vj 2 ^ is nearly identical to Since a large self-energy 

correction results in a short quasi-particle lifetime, 
suppresses T c while V { V enhances it, and ~ 3|U t (2) | 
in the Hubbard model generally results in a T c , if any, 
too low to be detected in FLEX calculations. [ po][27| ] This 
difficulty is resolved for large U c h and/or small Ujp , for 

which \Vtf > \ approaches V^. 

The microscopic origin of the charge fluctuation re- 
mains to be identified. Since there is no orbital degen- 
eracy in (TMTSF)2X, the origin cannot be the one pro- 
posed by Takimoto for Sr2Ru04. In fact, the mecha- 
nism of 2k F SDW-CDW coexistence in (TMTSF) 2 PF 6 
has been investigated by several authors. Some assume 
electron-lattice coupling, j2^j2^] while others envisage a 
purely electronic origin in terms of off-site repulsions up 
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